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Abstract 

We provide a quantum analysis of a DC SQUID mechanical displacement detector within the sub- 
critical Josephson current regime. A segment of the SQUID loop forms the mechanical resonator 
and motion of the latter is transduced inductively through changes in the flux threading the loop. 
Expressions are derived for the detector signal response and noise, which are used to evaluate the 
position and force detection sensitivity. We also investigate cooling of the mechanical resonator 
due to detector back reaction. 

PACS numbers: 85.25.Dq; 85.85.+j; 03.65. Ta 



1 



I. INTRODUCTION 



In a series of recent experiment s^^ and related theoretical work,-£>&L&>2 it was demon- 
strated that a displacement detector based on either a normal or superconducting single 
electronic transistor (SSET) can resolve the motion of a micron-scale mechanical resonator 
close to the quantum limit as set by Heisenberg's Uncertainty Principle.— The displace- 
ment transduction was achieved by capacitively coupling the gated mechanical resonator to 
the SSET metallic island. When the resonator is voltage biased, motion of the latter changes 
the island charging energy and hence the Cooper pair tunnel rates. The resulting modula- 
tion in the source-drain tunnel current through the SSET is then read out as a signature of 
the mechanical motion. 

Given the success of this capacitive-b&sed transduction method in approaching the quan- 
tum limit, it is natural to consider complementary, inductive-based transduction methods 
in which, for example, a superconducting quantum interference device (SQUID) is similarly 
used as an intermediate quantum-limited stage between the micron-scale mechanical res- 
onator and secondary amplification stages.— liiiiiLiS Unavoidable, fundamental noise sources 
and how they affect the SSET and SQUID devices are not necessarily the same. Furthermore, 
achievable coupling strengths between each type of device and a micron-scale mechanical 
resonator may be different. Therefore, it would be interesting to address the merits of the 
SQUID in comparison with the established SSET for approaching the quantum limit of 
displacement detection. 

In the present paper, we analyze a DC SQUID-based displacement detector. The SQUID 
is integrated with a mechanical resonator in the form of a doubly-clamped beam, shown 
schematically in Fig. [TJ Motion of the beam changes the magnetic flux $ threading the 
SQUID loop, hence modulating the current circulating the loop. We shall address the 
operation of the SQUID displacement detector in the regime for which the loop current is 
smaller than the Josephson junction critical current I c and at temperatures well below the 
superconducting critical temperature. We thus assume that resistive (normal) current flow 
through the junctions and accompanying current noise can be neglected. (See for exam ple 



Rcf. 



171 for a quantum noise analysis of resistively shunted Josephson junctions and Ref. 



18 



for a related analysis of the DC SQUID.) Such an assumption cannot be made with the usual 
mode of operation for the SSET devices, where the tunnel current unavoidably involves the 
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quasiparticle decay of Cooper pairs, resulting in shot noise. 

As noise source, we will consider the quantum electromagnetic fluctuations within the 
pump/probe feedline and also transmission line resonator that is connected to the SQUID. 
This noise is a consequence of the necessary dissipative coupling to the outside world and 
affects the mechanical signal output in two ways. First, the noise is added directly to the 
output in the probe line and, second, the noise acts back on the mechanical resonator via 
the SQUID, affecting the resonator's motion. 

With the Josephson junction plasma frequencies assumed to be much larger than the 
other resonant modes of relevance for the device, the SQUID can be modeled to a good 
approximation as an effective inductance that depends on the external current I entering 
and exiting the loop, as well as on the applied flux. In this first of two papers, we shall make 
the further approximation of neglecting the /-dependence of the SQUID effective inductance, 
which requires the condition I <C I c - I n the sequel,— we will relax this condition somewhat 
by including the next to leading 0(I 2 ) term in the inductance and address the consequences 
of this non-linear correction for quantum-limited displacement detection. 

Modeling the SQUID approximately as a passive inductance element, the transmission 
line resonator- mechanical resonator effective Hamiltonian is given by Eq. (|24p . This Hamil- 
tonian describes many other detector-oscillator systems that are modeled as two coupled 
harmonic oscillators, including the examples of an LC resonator capacitively coupled to a 
mechanical resonator— >2i and an optical cavity coupled to a mechanically compliant mirror 
via radiation pressure;— i2iL22»25a2£i the various systems are distinguished only by the depen- 
dences of the coupling strengths on the parameters particular to each system. Thus, many 
of the results of this paper are of more general relevance. 

The central results of the paper are Eqs. (1691) and (1701) . giving the detector response to a 
mechanical resonator undergoing quantum Brownian motion and also subject to a classical 
driving force. In the derivation of these expressions, we do not approximate the response 
as a perturbation series in the coupling between the SQUID and mechanical resonator as 
is conventionally done, but rather find it more natural to base our approximations instead 
on assumed weak coupling between the mechanical resonator and its external heat bath 
and weak classical driving force. Thus, in the context of the linear response paradigm, our 
detector should properly be viewed as including the mechanical resonator degrees of freedom 
as well, with the weak perturbative signal instead consisting of the heat bath force noise 
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and classical drive force acting on the mechanical resonator. Since the quality factors of 
actual, micron-scale mechanical resonators can be very large at sub-Kelvin temperatures 
(E.g., Q ~ 10 5 in the experiments of Refs. hid), quantum electromagnetic noise in the 
transmission line part of the detector can have strong back reaction effects on the motion of 
the mechanical resonator, even when the coupling between the resonator and the SQUID is 
very weak. One consequence that we shall consider is cooling of the mechanical resonator 
fundamental mode, which requires strong back reaction damping combined with low noise. 
Nevertheless, as we will also show, one can still analyze the quantum-limited detector linear 
response to the mechanical resonator's position signal using general expressions ( 1691) and 
(I70p ; under the appropriate conditions of small pump drive and weak coupling between the 
SQUID and mechanical resonator such that back reaction effects are small. 

The outline of the paper is as follows. In Sec. [Til we write down the SQUID-mechanical 
resonator equations of motion corresponding to the circuit scheme shown in Fig. [1] and then 
derive the Heisenberg equations for the various mode raising and lowering operators, subject 
to the above-mentioned approximations. In Sec. IHIl we solve the equations within the linear 
response approximation to derive the detector signal response and noise. In Sec. \1V\ we 
analyze both the position and force detection sensitivity, and address also back reaction 
cooling of the mechanical resonator. Sec. IVl provides concluding remarks. 



II. EQUATIONS OF MOTION 



A. Transmission line-SQUID-mechanical oscillator Hamiltonian 

Fig. [1] shows the displacement detector scheme. The device consists of a stripline resonator 
(transmission line T) made of two sections, each of length 1/2, connected via a DC SQUID 



(see Refs. 



27 



28 



29 



)30l for related, qubit detection schemes) . The transmission line inductance 
and capacitance per unit length are and Ct respectively. The Josephson junctions in each 
arm of the SQUID are assumed to have identical critical currents I c and capacitances Cj. A 
length Z osc segment of the SQUID loop is free to vibrate as a doubly-clamped bar resonator 
and the fundamental flexural mode of interest (in the plane of the loop) is treated as a 
harmonic oscillator with mass m, frequency uj m and displacement coordinate y. The total 
external magnetic flux applied perpendicular to the SQUID loop is given by <3> e xt + A-E> e xJosc2/, 
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where $ ex t is the flux corresponding to the case y = 0, B ext is the normal component 
of the magnetic field at the location of the vibrating loop segment (oscillator), and the 
dimensionless parameter A < 1 is a geometrical correction factor accounting for the non- 
uniform displacement of the doubly-clamped resonator in the fundamental flexural mode. 
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FIG. 1: Scheme for the displacement detector showing the pump/probe line l p\ transmission line 
resonator T', and DC SQUID with mechanically compliant loop segment having effective mass m 
and fundamental frequency to m . Note that the scale of the DC SQUID is exaggerated relative to 
that of the stripline for clarity. 



The transmission line is weakly coupled to a pump/probe feedline (p), with inductance 
and capacitance per unit length L p and C p respectively, employed for delivering the input and 
output RF signals; the coupling can be characterized by a transmission line mode amplitude 
damping rate 7 p t (see section III Bl below). Other possible damping mechanisms in the 
transmission line may be taken into account by adding a fictitious semi-infinite stripline 
environment (e), weakly coupled to the transmission line characterized by mode amplitude 
damping rate 7 e T-— While 7 e T can be made much smaller than 7 p t with suitable transmission 
line resonator design, we shall nevertheless include both sources of damping in our analysis so 
as to eventually be able to gauge their relative effects on the detector displacement sensitivity 
[see Eq. (1901)]. The SQUID, on the other hand, is assumed to be dissipationless. The 
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mechanical oscillator is also assumed to be coupled to an external heat bath (b), characterized 
by mode amplitude damping rate 7& m . 

A convenient choice of dynamical coordinates for the SQUID are j± = (0i ± fa) /2, where 
fa and fa are the gauge invariant phases across each of the two Josephson junctions.— For the 
transmission line, we similarly use its phase field coordinate fax, t),— where x describes 
the longitudinal location along the transmission line: — 1/2 < x < 1/2, with the SQUID 
located at x — 0. In terms of </>, the transmission line current and voltage are 

t ( +\ $ o dfax,t) 

lT ^ t) = -2^^x~ (1) 



and 



^OM)-^— ^— , (2) 



where $o = h/(2e) is the flux quantum. Neglecting for now the couplings to the feedline, 
stripline and mechanical oscillator environments, the equations of motion for the closed 
system comprising the superconducting transmission line-SQUID-mechanical oscillator are 
as follows (see, e.g., Ref. [14| for a derivation of related equations of motion for a mechanical 
rf-SQUID): 

^ " ( WS, (3) 



w, 2 7- + cos(7+) sin(7_) + 2f3 L l 



dt 2 dx 2 

($cxt + ^B cxt l osc y) 



7_ — 7T I n + 



o 



(4) 



^j 2 7+ + sin(7+) cos(7_) - ^ = 0. (5) 



L c 

and 

2 $o 



my + mu m y r A J B cxt / osc 7_ = 0, (6) 



where ujj = ^/2tiI c / (C j$o) is the plasma frequency of the SQUID Josephson junctions, 
the dimensionless parameter (3l = 27rL/ c /$o, L is the self inductance of the SQUID, n is 
an integer arising from the single-valuedness condition for the phase 27_ around the loop, 
and It is shorthand for It{x = 0, t). Eq. ([3]) is simply the wave equation for the phase field 
coordinate <f)(x, t) of the transmission line. Eq. (J3J) describes the current circulating the loop, 
which depends on the external flux threading the loop. Eq. (j3J) describes the average current 
threading the loop, which from current conservation is equal to one-half the transmission 
line current at x = 0. With the circulating SQUID current given by $ 7-/( 7r -^) (up to a 
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$ext dependent term), we recognize in Eq. ([6]) the Lorentz force acting on the mechanical 
oscillator. 

In addition to the equations of motion, we have the following current and voltage bound- 
ary conditions: 

I T (x = ±1/2, t) = (7) 

and 

where the external flux and current- dependent, effective inductance L e fr [$ ext (y), It] of the 
SQUID as 'seen' by the transmission line is 

^[*«t(y),/r] = ^g + j, (9) 

with $ext(?/) — ^ext + ^B e ^ t l osc y. Note that we have set n = 0, since observable quantities 
do not depend on n. 

We now make the following assumptions and consequent approximations: (a) ujj ^> lot 
u m (where lot is the relevant resonant mode of the transmission line); neglect the SQUID 
inertia terms toJ 2 j±. (b) (3l <C 1; solve for 7± as series expansions to first order in (3l- 
(c) \B ey± l osc y\ /$ 1; series expand the equations of motion to first order in y(t). (d) 
\I T /I C \ = 27 rLT-f 9 ^dx^ ^ ^> ser i es expand the equations of motion to second order in It- 

With cuj's typically in the tens of GHz, assumption (a) is reasonable. From Eq. (J4]), we 
see that a small (3l value prevents the 7_ coordinate from getting trapped in its various 
potential minima, causing unwanted hysteresis. With the 7+ expansion in It consisting of 
only odd powers, approximations (a) and (d) amount to describing the SQUID simply as a 
current independent, $ ext -tunable passive inductance element L eS [<3> ext (?/)] that also depends 
on the mechanical oscillator position coordinate y. Including the next-to-leading, I T term in 
the 7+ expansion gives an /^-dependent, nonlinear correction to the SQUID effective induc- 
tance. The consequences of including this nonlinear correction term for the quantum-limited 
displacement detection sensitivity will be considered in a forthcoming paper. 19 Solving for 
7 + to order It and substituting in Eq. ([9]), we obtain: 

L cS [$ext(2/)] « ^ sec [—^- ) , (10) 

where the self inductance L contribution has been neglected since it is of order (3l <C 1. 
Solving for 7_ to order I T and substituting into Eq. (jSJ), we obtain for the mechanical 



oscillator equation of motion: 

my + muj 2 m y - 7rAB ° rtZ ° ac/ r tan ^ cxt /^ Q ) sec (7r$ cxt /$ ) = 0, (11) 



where from (c), we have set y = in the solution for 7_ and have dropped an overall 
constant term. Since the 7_ expansion in It consists only of even powers, we must go to 
second order in It so as to have a non-trivial transmission line-oscillator effective coupling. 
Thus, the SQUID phase coordinates j± have been completely eliminated from the equations 
of motion, a consequence of approximation (a); the SQUID mediates the interaction between 
the mechanical oscillator coordinate y and transmission line coordinate without retardation 
effects. 

From Eq. (TIT]) , it might appear that the force on the mechanical oscillator due to the 
transmission line can be made arbitrarily large by tuning $ ext close to $o/2. Note, however, 
that the proper conditions for the validity of the It and (3l expansions are: 

I T 



SeC (7T$ ext /$o) 
J-r. 



< 1 (12) 



and 

|/3 L sec(7r$ cxt /$o)| < 1. (13) 

We now restrict ourselves to a single transmission line mode and derive approximate 
equations of motion for the mode amplitude. Suppose that the mechanical oscillator position 
coordinate is held fixed at y = 0. The following phase field satisfies the current boundary 
conditions (JTj): 

f -<f>{t)co8[ko(x + l/2)]; x<0 
<f>(x,t) = < , (14) 

I +(f>(t) cos [ko {x - 1/2)} ; x > 
with the wavenumber fc determined by the voltage boundary condition (jSJ): 

kf)l ( kol\ LtI 



tan ^ = - - ;* - . (15) 
2 V 2 J L eS ($ ext ) 1 ; 



The wave equation ([3]) gives for the transmission mode frequency: lot = ko/ \JLtCt- Substi- 
tuting the phase field ( |T4T) into the It part of the oscillator equation of motion ( TTTi) further- 
more gives the transmission line force acting on the oscillator with fixed coordinate y = 0. 
Now release the mechanical oscillator coordinate and suppose that for small [condition (c)] , 
slow [condition (a)] displacements, the force is the same to a good approximation. Then the 
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oscillator equation of motion becomes 



my(t) + mw 2 m y{t) + X -C T l ) sin 2 (k l/2) 



x 







tan (7r$ext/$o) sec (7r$ ex t/$o) 



u%<ff(t) = 0, 



(16) 



($ /2tt) 4ttL t // c 

From Eq. ( fl6l) . we can determine the mechanical sector of the Lagrangian, along with the 
interaction potential involving y and the mode amplitude 6. The remaining transmission 
line sector follows from the wave equation ([3]) and we thus have for the total Lagrangian: 

/ . > 1 „ 1 „ „ 1 /$ n \ 2 
L \(J),y,4>,y , 



x <^ -dr 



\mi? - \muly 1 + l -C T l (^J sin 2 (k l/2) 
1 - ■ t-?^- tan (7T$ cxt /$ ) sec (7r$ cxt /$ ) 



2 J.2 
UJ T (p 



(17) 



($ /2tt) 4ttL t Z/ c 

From Eq. (fl7|) . we see that for motion occuring on the much longer timescale u;" 1 ^> oo^, 
the mechanical oscillator has the effect of modulating the frequency of the transmission line 
mode. 

The associated Hamiltonian is 

2 



C T Z(f*) 2 sin 2 (M/2) 

■ T~T~TT tan (**W*o) sec (vr$ cxt /$ 
($o/27r) 4itL t II c 



" 2 j 2 
-UJ T (p 



Py 2 2 

2m 2 mU 



(18) 



Let us now quantize. For the transmission line mode coordinate, the raising(lowering) 
operator is defined as: 

2 



1 



'2Huj t [\C T l ($o/2tt) 2 sin 2 (k l/2)] 
and for the mechanical oscillator 

1 



1 / $ 

2^ 1 2~7r ^ ^°^ //2 ' ) Wt< ^ ^ 



(19) 



(muy =F zp y ) . 



(20) 



y2mujh 

In terms of these operators, the Hamiltionian (|18p becomes (for notational convenience we 
omit from now on the 'hats' on the operators and also the 'minus' superscript on the lowering 
operator): 

1 2 

H = hu) T a\a T + hu) m a%fl m + -huj T K Tm (a T + of) (a m + a+) , (21) 



where the dimensionless coupling parameter between the mechanical oscillator and trans- 
mission line mode is 

K Tm = - XB ^7^ ZP A ru tan (*W*o) sec (7r$ cxt /$ ) , (22) 

with Ax zp = y/hf (2mw m ) the zero-point uncertainty of the mechanical oscillator. From 
expression (jTOl) for the effective inductance, another way to express the coupling parameter 
is as follows: 

XB ext l osc Ax zp $ dL eS /d§ 

ext /no\ 

($ /27t) 7T L T / 

From Eq. ( |23|) . we see that in order to increase the coupling between the mechanical oscillator 
and transmission line, the SQUID effective inductance-to-transmission line inductance ratio 
must be increased. The advantage of using a SQUID over an ordinary, geometrical mutual 
inductance between a transmission line and micron-sized mechanical oscillator is that the 
former can give a much larger effective inductance. As we shall see in Sec. IIV} just requiring 
that the inductances be matched such that go dL cS /d^ c ^t ^ \ [ s sufficient for strong back 
reaction effects with modest drive powers, even though the other term in Kttu describing 
the flux induced for a zero-point displacement is typically very small. 

Assuming then that Kj; m -C 1 and making the rotating wave approximation (RWA) for 
the 'T' part of the interaction term in the system Hamiltonian (12T1) . i.e., neglecting the terms 
(ax) 2 and (a^) 2 , we have (up to an unimportant additive constant): 

H = hjj T a^a T + ^a; m a+a m + hjj T K Tm a^a T (a m + a+) . (24) 

Many other systems are modeled by this form of Hamiltonian, a notable example being 
the single mode of an optical cavity interacting via radiation pressure with a mechanically 
compliant mirror.— Thus, much of the subsequent analysis will be relevant to a 
broad class of coupled resonator devices-not to just the transmission line- SQUID- mechanical 
resonator system. 



B. Open system Heisenberg equations of motion 

So far, we have treated the transmission line and mechanical resonator as a closed system 
with SQUID-induced effective coupling . Of course, a real transmission line mode will expe- 
rience damping and accompanying fluctuations, not least because it must be coupled to the 
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outside world in order for its state to be measured. Furthermore, the mechanical resonator 
mode will of course be damped even when decoupled from the SQUID. It is straightforward 
to incorporate the various baths and pump/probe feedline in terms of raising/lowering op- 
erators. Assuming weak system-bath couplings, which again justify the RWA, we have for 
the full Hamiltonian: 

H = hwra^aT + huo m a^ n a m + hiOTK Tm a^,aT (a m + a+) 

+h J duouoa p (co)a p (co) + h J duouoap {uo)a e {uo) + h J duouoap (uo)a b (uo) 

+h J duo \K* T a^(uo)a T + K pT a^a p (uo)] + h J duo [K* T a^(co)aT + K pT a^a e (co)] 
+h j duo [K bm a+(uo)a m + K bm a+a b (uo)] - \J 2 ^ («m + a+)F cxt (t), (25) 

where a p denotes the pump/probe (p) feed line operator, a e the transmission line bath ('e' for 
'environment') operator, and a b the mechanical resonator bath (6) operator. These operators 
satisfy the usual canonical commutation relations: 

[oi(u),at(u')] =5 t ,5(uo-uo'). (26) 

The couplings between these baths and the transmission line and mechanical resonator 
systems are denoted as K p t, K b t, and K bm . Note we have also included for generality 
a classical driving force F ext (t) acting on the mechanical resonator. This allows us the 
opportunity to later on analyze quantum limits on force detection in addition to displacement 
detection. 

Within the RWA, it is straightforward to solve the Heisenberg equations for the bath 
operators and substitute these solutions into the Heisenberg equations for the transmission 
line and mechanical oscillator to give 



da m . i h + 

—— = -iuo m a m + - 4 /- F ext (t) - iuo T K Tm a T a T 

dt n\ 2muo m 

-J du\K Tm \ 2 J dt'e'^-^a^-i J duoK bm e-^-^ 



a b (uo,t ) (27) 



and 



dar . T , / +\ 

—j- = -iuo T a T - iuo T K Tm a T [a m + a m ) 

- [ duo\K pT \ 2 [ dt'e-^WarW-i f duoK pT e- tu)(t - to) a p (uo,t ) 

J Jt J 
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- J du \K eT \ 2 J dt'e-^-^ciTit') - i j duK eT e- iu}{t ~ to) a e {uj, t ). (28) 

We now make the so-called 'first Markov approximation',—^ in which the frequency 
dependences of the couplings to the baths are neglected: 



K pT {u) = J^fe^ 



K eT {uj) = \\— e 



7eT _j0 eT 



KU^) = J 7 ^***", (29) 

where the 7's and 0's are independent of u as stated. The Heisenberg equations of motion 
(1271) and (j28) then simplify to 



= ~iu m a m + \\ 7—^ — F ext (t) - iuj T K Tm a^a T 
at ft V 2mu; m 

-% m a m (t) - i^Ze^atit) (30) 



and 



ddT 
~dT 



-iu T a T - iu T K Tm a T (a m + a+) 
- lpT a T {t) - iy/^e^^t) 



-le T a T (t) - iy/2^+«-c£(t), (31) 

where the 7i's are the various mode amplitude damping rates (assumed much smaller than 
their associated mode frequencies) and the 'in' operators^^ 1 ^ 1 ^ are defined as 

aj n (t) = -L= [ due'^-^aiiu^o), (32) 



/2tt 

with t > to. The time to can be taken to be an instant in the distant past before the 
measurement commences and when the initial conditions are specified (see below). We can 
similarly define 'out' operators: 

a° ut (t) = — L / dwe-W^aiiuth), (33) 
v2vr J 

with t\ > t. The time t\ can be taken to be an instant in the distant future after the measure- 
ment has finished. From the Heisenberg equations for the bath operators and the definitions 
of the 'in' and 'out' operators, we obtain the following identities between them:— ^ 



a° ut (t)-<(t) = -zV^e- l ^a T (t) 
12 



< Ut (t)-a[ n (t) = -iy/^e-^Ornit) 

a^{t)-af{t) = -iy/^e-^arit). (34) 

In outline, the method of solution runs in principle as follows:— iMJ&M. (i) specify the 'in' 
operators. (2) Solve for the system operators a m (t) and ar(i) in terms of the 'in' operators. 
(3) Use the relevant identity fl34|) to determine the 'out' operator a° nt (t), which yields the 
desired probe signal. It is more convenient to solve the Heisenberg equations in the frequency 
domain with the Fourier transformed operators 0{t) = ^= doje~ lult O{uj). The equations 
for the system operators then become 

1 



a m (u) = — { y/2 lbm e^af(iu) - F ext (u) 

+ uJ I Krn L f°° dJ ^jy+^j + a+^ariu + u')] \ (35) 



2V2tt J- 

and 



a T (co) = -± Jy^e^a^ + y^e^^iu) 



OC 



Tm I du'a T (uj') [a m (uj - u') + a+(o>' - u)] [• , (36) 



/2tt 

while the relevant 'in/out' operator identity becomes 



•y^e-^wriu) + aj». (37) 



C. Observables and 'in' states 

Before proceeding with the solution to Eqs. (135j) and (1361) . let us first devote some time 
to deriving expressions for observables that we actually measure in terms of a° ut (u;). Model 
the pump/probe feedline as a semi-infinite transmission line — oo < x < 0. Solving the wave 
equation for the decoupled transmission line and then using the expressions (pQ), ([2]) relating 
the current /volt age to the phase coordinate, we obtain 




and 



I™\x,t) = - / duJ-—mn(ux/v p ) [e~ iwt 0°^ (u) + e^a° ut+ (^)] (38) 



/OO / ry p- 

du y =^ cos (ujx/v p ) [e~ iult a° ut (w) - e^'a° ut+ (a;)] , (39) 



13 



where the sinusoidal x dependence in the current expression follows from the vanishing of the 
current boundary condition at x — 0, the feedline impedance is Z p = L p /C p and the wave 
propagation velocity is v p = 1/ ^JL V C V . Suppose the current/volt meter is at x — > — oo, so 
that the actual observables correspond to measuring the left-propagating component of the 
current/voltage. Then decomposing the x-dependent trig terms into their real and imaginary 
parts, we can identify the left propagating current /voltage operators as 



I out (x,t) = 

and 

V out (x,t) 




The output signal of interest due to the mechanical oscillator signal input will lie within 
some bandwidth Soj centered at u s , the 'signal' frequency, and so we define the filtered 
output current J out (x,t\uj s ,Su) and voltage V ont (x,t\uj s ,Su) to be the same as the above, 
left-moving operators, but with the integration range instead restricted to the interval 
[uj s — 5u/2, uj s + 5u/2]. 

Since the motion of the mechanical resonator modulates the transmission line frequency, 
one way to transduce displacements is to measure the relative phase shift between the 'in' 
pump current and 'out' probe current using the homodyne detection procedure.— Another 
common way is to measure the 'out' power relative to the 'in' power, or equivalently the 
mean-squared current /volt age (all three quantities differ by trivial factors of Z p ). We will 
discuss the latter method of transduction; the former, homodyne method can be straight- 
forwardly addressed using similar techniques to those presented here. Thus, we consider the 
following expectation value: 

( [SI out (x, t\u a , Su)] 2 ) = ( [l out (x, t\u s , 5u)} 2 ) - (l out (x, t\u„ 5uj)) 2 , (42) 

where the angle brackets denote an ensemble average with respect to the 'in' states of 
the various baths and feedline (see below). If the mechanical oscillator is being driven 
by a classical external force whose fluctuations are invariant under time translations, i.e., 
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{F ext (t)F eyLt (t')) = C(t — t'), then the above, mean-squared current will be time- independent. 
Alternatively, if F ext (t) is, e.g., some deterministic, AC drive, then we must also time-average 
so as to get a time-independent measure of the detector response: 

([8I^(x,t\u s ,8uj)} 2 ) = — / dt([l 0Ut {x,t\u s ,Su)] 2 ), (43) 

J-mJ-t m /2 x 1 

where Tm is duration of the measurement, assumed much larger than all other timescales 
associated with the detector dynamics. We have also assumed that the time-averaged cur- 



rent vanishes in the signal bandwidth of interest: (J out (u s , 8oj)) = 0. Substituting in the 
expression (HP!) for J out (x, t\u s , 80S) in terms of the a° ut operators, we obtain after some 
algebra: 

( Wa ,M] 2 > = ±r r^^^h.A- 2 — -sin [fa- u 2 )T M /2]\ 

x i (d^Mc^M + a^Wa%*M) ■ (44) 

As 'in' states, we suppose ksT -C fkox, such that the relevant transmission line 'in' bath 
modes (u e ~ ujt) are assumed to be approximately in the vacuum state. On the other hand, 
with the mechanical mode typically at a much lower frequency u m <C lot, we assume that its 
relevant 'in' bath modes (oof, ~ u m ) are in the proper, non-zero temperature thermal state. 
For the pump/probe feedline, we consider the following coherent state: 30 



|{a(w)}) p = exp 
where |0) is the vacuum state and 



dua(u;) (af + (uj) -c#(u})) 



|0L, (45) 



^"-w^ 5 — • (46) 



normalized such that the amplitude of the expectation value of J in [the right propagating 
version of (HOI) with a° ut replaced by a™] with respect to this state is just I . Again, we 
suppose kftT frw p , so that thermal fluctuations of the feedline are neglected. The fre- 
quency width of this pump drive is assumed to be the inverse lifetime of the measurement. 
Below we shall see that the output mechanical signal will appear as two 'satellite' peaks on 
either side of the central peak at u p due to the pump signal, i.e, the mechanical signal can 
be extracted by centering the filter at either of uo s = u p ± u m (up to a renormalization of 
the mechanical oscillator frequency), corresponding to the anti-Stokes and Stokes bands. 
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Note that we do not have to specify the initial to states of the mechanical resonator and 
transmission line systems; ay(io) and a m (to)-dependent initial transients have been dropped 
in the above equations for a<r(w) and a m (uj), since they give a negligible contribution to the 
long-time, steady-state behavior of interest. 



III. SOLVING THE EQUATIONS OF MOTION 



A. Linear response approximation 



We are now ready to solve for <([5/ out ] Introduce the following shorthand notation: 
S t (uj) = v^e^aJV) + v^e^4 n H 



K 



ujtK> 



Tm 



(47) 



and 7r = 7 P T + 7eT, the net transmission line mode amplitude dissipation rate due to loss via 
the probe line and the transmission line bath. Substituting Eq. (1351) for a m (uj) into Eq. (|36|) 
for ax(uj) yields the following, single equation in terms of ax(uj) only: 



a T (u) = / duj'a T (uj — uj')A(u, uj') + / duj'B(uj, uj')o,t{uj — uj') 

J — oo J —oo 

/oo 
duj" [a T (uj")a+(uj" - uj') + a+(uj")a T (uj" + uj')] + C(u), (48) 
-oo 

where, for the convenience of subsequent calculations, we have made this equation as concise 
as possible with the following definitions: 



A(uj, uj') 
B(uj,uj') 
C(oj) 



UJ — LOt + «7t 

/C 2 /2 

uj — ujt + i^T 
St(uj) 



S m (uj') 



UJ' - UJ m + i-fbr, 



+ 



-uj 



-UJ — UJ r , 



1 



uj' - uj m + i^ bri 



-UJ — UJ r , 



- -■ (49) 

uj — uj t + 

We expand Eq. fj48|) for a^(^) to first order in the mechanical oscillator bath operator 
a^(uj) and external driving force F ext (uj) [equivalently expand in A(uj, uj')}: cit{uj) ~ a^\uj) + 
a^\uj), where 



aip (uj) 



duj'B(uj,uj')a^\uj — uj') 
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doo" 



aP (u")aP + (uo" - uj') + 4 0)+ (o/')4 0) (a/' + oo')} + C(w) (50) 



and 



a^'(u;j 



x 



dw / a^ ) (u;-w / )^(w,a/)+ / du/B(u,<jf)c$\u - oo') 

30 

da/' 



— oo 

oo 



(0)/ //x (0)+/ // i\ . (0)+/ //x (0)/ // . /x 

a T (uo )a T ' [uo — uj J + a T [oo )a T [oo + uo ) 

/oo 
da/' cip (uj")a!p + (uj" — uj') 
■oo 

. (1)+/ //\ (0)/ // . a . (1)/ //x (0)+/ // /x 

(a; Ja^ (a; + lo J + a T (a? Ja T (c<j — c<j J 

, (0)4- / //\ (1)/ // , a" 

+a T (uj )a T [uj + uj J 



(51) 



Eq. (1501) then yields the detector noise, while ( 15"TI) yields the detector response to the signal 
within the linear response approximation. Thus, our approach here is to treat the mechanical 
oscillator as part of the detector degrees of freedom, with the signal defined as the thermal 
bath fluctuations and classical external force acting on the oscillator. This is the appropriate 
viewpoint for force detection. On the other hand, if the focus is on measuring the quantum 
state of the mechanical oscillator itself, then the oscillator should not be included as part 
of the detector degrees of freedom. Nevertheless, as we shall later see, the latter viewpoint 
can be straightforwardly extracted from the former under not too strong coupling K^ m and 
pump drive current amplitude Iq conditions. 



B. Semiclassical approximation 

The sequence of solution steps to Eqs. (150]) and (l5Tj) are in principle as follows: (1) Solve 
first equation (150|) for aP(uj) in terms of B(u,uj') and C(uj); (2) Substitute the solution for 
aP {uj) into Eq. (l5"Ti) for a^ (oo) and invert this Eq. (which is linear in a^" (oo)) to obtain 
the solution for (oo) in terms of A(uj,uj'), B(uj,uj'), and C(oo). It is not clear how to 
carry out these steps in practice, however, since the equations involve products of non- 
commuting operators. Thus, we must find some way to solve by further approximation. 
The key observation is that the feedline is in a coherent state, which is classical-like for 
sufficiently large current amplitude I so as to ensure signal amplification. We therefore 
decompose (oo) into a classical, expectation- valued part and quantum, operator- valued 
fluctuation part, aP(oo) = /aP(uj)\ + 5aP(uo), and subtitute into Eq. ( 1501 for af/ (oo), 
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linearizing with respect to the quantum fluctuation ddj, (oo). This gives two equations, one 
for the expectation value 



doo'B(oo,oo') {0^(00-00') 



duo" 



a^' [oo ))\cl t {oo — uo ) 



+ (4 0) V) > < 4°V +«>')) + (CM) 



aP (oo) 



and the other for the quantum fluctuation: 

5aP(uo) = j doo'B(u,u')8af(oo - uo') j doo" ^4°V')) (4° )+ M' 
+ ' 

x 



(52) 



oo 



4° )+ M')) (4°V + "'))] + j duo'B(uo,uo') (a { °\uo - uo')} 
5af(uo") (aP + {uo" - uo')} + (4°V')} <*4° )+ M' - 



+<54° )+ (c/') (af(oo" +<J)) + (4° )+ M')} 54 0) M' + + <*C(w). (53) 

Eq. (EH) for 4^M i s approximated by replacing a^\oo) with its expectation value (a^\oo)}, 
i.e., we drop the quantum fluctuation part 8a!j)(oo). This is because Eq. ( IS"Tp already depends 
linearly on the quantum fluctuating signal term A(uo,uo'), which we of course want to keep. 
Dropping the 5a^\uo) contribution to Eq. fl5Tj) amounts to neglecting multiplicative detector 
noise, which is reasonable given that we are concerned with large signal amplification. 



C. Complete solution to detector signal response and noise 

The sequence of solutions steps are therefore in practice as follows: (1) Solve Eq. (1521) first 
for (^a^\uo)}; (2) Substitute this solution into Eq. (|5"TT) for a^\uo) and invert; (3) Substitute 
the solution for (^a^ (uo)} into the Eq. (1531) for 5a^\uo) and invert; (4) Use these solutions for 
a^\uo) and 5a^\uo) to determine the detector signal and noise terms, respectively. Beginning 
with step (1), we have 



i 



t2 



(C(oo)) = - v m \ «(u,)) = v m \ ■ lJ^L e -^ lT Ml\ (54) 

where Auo = uo p — uot is the detuning frequency (not to be confused with the bandwidth 5uo) 
and note (a^u;)) = (recall, we assume the transmission line resonant frequency oot mode 
is in the vacuum state). Given that Tm is the longest timescale in the system dynamics, 
(C(oo)) is sharply peaked about the frequency uo p and we will therefore approximate the 
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exponential with a delta function: (C(uj)) = c5{uj — uj p ), where 



27re*w \IlZ plpT ^ 



7^ — iAuj y hjjj. 



v 



Considering for the moment an iterative solution to Eq. (|52p for (a^\uj)), we see that 
(<4? (oj)) must also have the form of a delta function peaked at uj p : (tZj? (a>)) = X'H^ — ^p)- 
Substituting this ansatz into Eq. (1521) . we obtain the following equation for x' 



X = 2 X \x\ B(u p ,0)+c (56) 

This equation has a rather involved analytical solution. For sufficiently large |c| 2 \B(u p , 0)| 
the response can become bistable (i.e., two locally stable solutions for x). This region will 
not be discussed in the present paper, however. When we consider actual device parameters 
later in Sec. IIVI we will assume sufficiently small drive such that X ~ c i allowing much 
simpler analytical expressions to be written down for the detector response. 

Proceeding now to step (2), we substitute the expectation value (^a^\uj)^ = x^( u ~ ^p) 
for the operator a^\uj) into Eq. ([5"T| for a^(u). Carrying out the integrals, we obtain 

{1-2 |x| 2 [B(w,0) + B(w,w- u p )]} a^\uj) - 2 X 2 B{uj,uj - uj p )a^ )+ (2uj p - uj) 

= xA(u,u-u p ). (57) 

Before we can invert to obtain a^\u), we require a second linearly independent equation 
also involving a^ + (2uj p — uj) and cSp (uj). This equation can be obtained by replacing u> with 
2uj p — uj in Eq. (157)1 and then taking the adjoint: 

{1 + 2 |x| 2 [B(u - 2Auj, 0) + B(uj - 2Auj, uj — uj p )]} a^ + (2uj p — uj) 
+2 X * 2 B(uj-2Auj,uj-uj p )cSP{uj) = ~ X *M^ ~ 2Auj, uj - uo p ), (58) 

where we have used the identities A + (2uj p — uj, uj p — uj) = — A(uj — Auj,uj — u p ), B*(2uj p — 
uj, uj p — uj) = —B{uj — 2Auj, uj — uj p ), and B*{2uj p — uj, 0) = —B(lj — 2Auj, 0). Inverting, we 
obtain 

(uj) = ai(uj)A(uj, uj — uj p ) + a2(uj)A(uj — 2Auj, uj — uj p ), (59) 

where 

ai(u) = Diu)- 1 ^ + 2 \ X \ 2 [B(uj - 2Auj, 0) + B(u - 2Auj,uj - uj p )}} X (60) 
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and 



with determinant 



a 2 (uj) = -2D(u) 1 \x\ 2 B(uj,uj - lu p )x 



(61) 



D(u) = {1 - 2\ X \ 2 [B(w,0) + B(uj,uj - uj p )}} 



x 



{ 1 + 2 |x| 2 [B(u - 2Auj, 0) + B(uj - 2Au, uj - uj p )} } 



+4 |x| 4 B(u, uj — uj p )B(uj — 2Au, uj — uj. 



pi- 



(62) 



Moving on now to step (3), we substitute the expectation value (o^ (u)) = X&i 00 ~~ ^p) 
into Eq. (153]) for 8a§ (o>) and carry out the integrals to obtain: 

{1-2 |x| 2 [B(w,0) + B(w,u- Up)}} 5aP(u) - 2 X 2 B(uj,uj- uj p )5a^ )+ (2uj p - uj) 

= 5C{uj). (63) 

Replacing uj with 2uj p — uj in Eq. (|63|) and then taking the adjoint: 



{1 + 2 |x| 2 [B(u - 2Auj, 0) + B(uj - 2Au, uj - uj p )} } 5a§ )+ (2u p - uj) 
+2 X * 2 B(uj - 2Auj, uj - uj p )5af(uj) = SC + (2uj p - uj). 

Inverting Eqs. fl63|) and fl64l) . we obtain 

5aP{uj) = f3 1 {uj)5C(uj) + f3 2 {uj)5C + {2uj p - uj), 



(64) 



(65) 



where 



and 



ft (w) = Diuj)- 1 ^ + 2 |x| 2 [B(uj - 2Auj, 0) + B(uj - 2Auj,uj- uj p )]} 



(66) 



(3 2 (uj) = 2D(ujy 1 X 2 B(uj, uj - uj p ) (67) 

We are now ready to carry out step (4). To obtain the detector response, we substi- 
tute into expression (j44p for ^[<5/ out ] 2 ) the linear response approximation to the 'out' probe 
operator [see Eq. (I37|) ]: 



-i^j2^ pT e 



-i<t> v T n ( l ) 



drr [UJ \ 



+ 



-i^fe^e-^SaPiu) + 6a p n (uj) 



(68) 



The first square-bracketed term will give the signal contribution to the detector response, 
while the second bracketed term gives the noise contribution. Note that the average values 
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,(°), 



LO 



and (a™ (a;)) are not required in the noise term since they give negligible contri- 



bution in the signal bandwidths of interest centered at lo s 
signal part of a° nt (u>), we obtain after some algebra: 



u p ± lo % 



Substituting in the 



([51°°* (oo s ,5uo)] 2 ) 
It 

OlAlO) 0(2(00 

X 



signal 
IpT 



u s +6u,/2 duJ r u 



1 T 



) 7 f + AlO 2 J Us _ 5u} /2 27T 

to — lo p + Aio + i-f T 



c 



c 

2^bm 



Up (tO — LOp + Alo) 2 + 7y 
2 



a; 



Mr, 



10, 



lo — LOp — Alo + ijT 



[2n(u - ujp) + 1] + 



bin 



2 + il, 



+ 



loKx-m^T 

It 



LO - LO m ) 2 + H m 



\OJr, 



[2n(co p — lo) + 1] 



7y + Alo 2 2mhio m ^jb 



a\\LO) a 2 (io) 



, sin [(a; - lo')T m /2] 
(lo - lo') T M /2 

276m 



lo — LOp + Alo + 
lo — LOp — Alo + i 7 T 

27fe m 



m Juj s ~Slo/2 
2 



LO 



It 



2n \_LOp (lo — LOp + Alo) 2 + 7 |> 



(LO -LOp- LO m ) 2 + 7^ 
-F e xt(^p - U;)F £ 



ext 



.to 



LOr, 



, . :tA a;p-a/j>. (69) 

(a;,, - u - Lo m y + 7bm j 

where n(a;) = (e hu) ^ kBT — l) 1 is the Bose-Einstein thermal occupation number average for 
bath mode lo. The signal part of the detector response comprises a thermal component 
and a classical force component. In the limit of weak coupling K^m — ► and or small 
drive current amplitude I — > 0, we have a.\(uj)/c — > 1, a 2 (Lo)/c — > and we note that the 
frequency resolved detector response has the form of two Lorentzians centered at lo p ± to m . 
The resulting expression for the detector response coincides with an 0(K^ m ) perturbative 
solution to the detector response (jHJ) via the linear response Eqs. ( 1501 and (l5Tj) (but no 
semiclassical approximation). However, as shall be described in Sec. IIV1 when the current 
drive is not small and or coupling is not weak, then the a$ terms will modify this simple 
form, at the next level of approximation renormalizing the Lorentzians, i.e., shifting their 
location and changing their width. 

Substituting in the noise part of a° nt (Lo), we obtain after some algebra: 

27t7pT 



out 



LO, 



MY) 



— rko- — 

u, s -Su,/2 2tt (lo - LOp + Alo) 2 + ^ 



x { IAMI'+ ^- tt> + ^! + 1 lAMI'-if|ftM] + ( ""^ + A " ) to[ft M ] 



[LO 



+z, 



_ x fko s 5lo 
~2~2k' 



LOp - Alo) 2 + 7 f, 



7t 



(70) 
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The noise part of the detector response comprises a back reaction component (the integral 
term) where transmission line noise drives the mechanical oscillator via the SQUID coupling, 
and a component that is added at the output due to zero-point fluctuations in the probe 
line. While not as obvious given the form of Eq. (|70|) . one may again verify (see Sec. lIVp that 
the detector back reaction on the mechanical oscillator takes the form of two Lorentzians 
centered at u p ± u m in the weak coupling and or weak current drive limit, coinciding with 
an 0(K^ m ) perturbative calculation. 

Eqs. ( 1691) and (ITU]) are the main results of the paper, their sum giving the net output 
mean-squared current. 



D. Quantum bound on noise 



As articulated by Caves,— the fact that the 'in' and 'out' operators satisfy canonical 
commutation relations places a lower, quantum limit on the noise contribution to the detector 
response, Eq. (1701) . We now derive this quantum limit. First write the 'out' operator (168|) 
as 

a° u » = -iy^e'^aPiu) + N(id), (71) 

where N(id) = -i^q^e'^? 'dap \u) + 5a™ (uj) is the noise part. Taking commutators, we 
have the following identity relating the noise and signal operator terms: 



[N(u), N + (cu')] = 5(uj - u') - 2 lpT \a^(id), a? + (cu' 



(72) 



Now, from the Heisenberg Uncertainty Principle, one can derive the following general in- 
equality: 

(N[f]N + [f]+N + [f]N[f]) > \{[N[f),N + [f)])\, (73) 

where N[f] = J °° duf(uj)N(uj) and f(uj) is an arbitrary function. Inserting the commutator 
identity (1721) . Eq. ( 1731) becomes 



(N[f]N + [f]+N + [f]N[f})> 



^|/H| 2 -2 7pT ([aS! ) [/],ai! )+ [/] 



(74) 



Choosing the 'filter' function f(u) = u)Q(oj — oj s + 5u)/2)Q(oj s + 5oj/2 — u) and evaluating 
the commutator, we obtain the following lower bound on the detector noise: 



([61°* (lo s ,6u)} 2 ) 



> 
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_ x HuJ s 5uJ 

IT 2^ 



I K Tm uj T \ 2 7pr f 
7T ; 7r + ^ 2 J u 



w .+&,/2 duJ 



c 



uj — u p + Auj + ijT 



LU — UJ 



J p 

p - Auj + i"Y T 

^^bm 



UJ S — 8uj 1 2 
2 



2tt 



UJ 



1 2 T 



oj p (uj — Up + Auj) 2 + 7|, 



(75) 



_ (u - uj p - uj m ) 2 + 7 fc 2 m (wp - w - uj m ) 2 + 7 fc 2 „ 

In the next section we will address the extent to which the detector noise can approach the 
quantum bound on the right hand side of Eq. (1751) . depending on the current drive amplitude 
Jo and other detector parameters. 



IV. RESULTS 



A. Analytical approximations 



To gain a better understanding of the detector response, we now provide analytical ap- 
proximations to Eqs. (I69I) and (170]) that are valid under the condition |c| 2 \B(u p , 0)| <C 1 
such that x ~ c [ see Eq. (|56|) ]. i.e., the expectation value (a^( 



v. 



for the transmission line 



depends approximately only on the pump/probe feedline state and not on the mechanical 
oscillator state. Explicitly, this condition reads: 

2IlZ v K 2 Tm uj Tl p T 



<1, 



(76) 



hu m (7! + Au 2 ) 3/2 

placing an upper limit on I and ifrm for the validity of this approximation. We also 
assume that the mechanical and transmission line mode frequencies are widely separated: 
u m <C ut, and with small damping rates: ■jbm <C w m i 7t <C wt> We do not restrict 
the relative magnitudes of uj m and 77-, however. A simple picture emerges in which the 
detector back reaction 'renormalizes' the mechanical oscillator frequency and damping rate: 
u m -> Riu^m and % m -> -R 7 7 bm , where 



u) m + ( Auj + 



C LUrpJ\rp 



Tm 



TTU) r , 



c\ 2 u^K 2 m [ri + Au 2 



ujz. 



7T [7| + (Aw + UJ m ) ] [jr + ( AuJ ~ Urn) ] 



(77) 



and 



n 



Ibra - \ Auj + 



7YUJ ri 



2\c\ 2 UJ%K% m uj m y T 
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7T [7I + (Auj + uj m ) ] [y 2 + (Auj - u m ) ] 



(78) 
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where c is defined in Eq. (1551) . With the measurement filter bandwidth centered at either of 
oj s = uj p ± RujUm, the approximation to Eq. (169]) for the signal response is (with the classical 
force term omitted): 

' IoK Tm u T \ 2 IpT 7t 



It J Tt + Auj2 1 T + ( AuJ ± 



([<5/ out (u s = lu p ± R„U} m , 5uj)} 2 ) = I 

signal V 

x / ^ 7 R 2Mm ^ + 1] • (79) 

When there is a classical force acting on the mechanical oscillator, we must add to Eq. ( 179]) 



the term 

2 ' ' "^ +5 ^ 2 duoduo' 7 2 



{ IoK Tm (jJ T \ IpT 

V It ) 7f + A 



Au; 2 2mh J j m rtbm J w a -6<v/2 2?r (^ - uj p + Acj) + 7f, 
sin [(w - u/)T A/ /2] J 2 7bm „ / \p* / / \ 

(W - £«/) T M /2 \ (a, - U„ - i^ m ) 2 + extV p, extl W 

^(wp-w^^-wOi. (80) 



(a> p - a; - R^m) 2 + (-R 7 7bm) 
The approximation to Eq. (I7D|) for the detector noise is 

'IoK Tm uj T \ 2 jpT 7f 



7T y 7y + Auj 2 ^2 + (Acj ± cu m ) 

^ S+ " /2 ^ *»= ,iV ± + Z^^, (81) 

.,-4w/2 2vr ( w - uj p =f i?^ m ) 2 + (i? 7 7 bm ) 2 2tt 

where the back reaction noise parameter is 

AT _ \c\ 2 K 2 m u 2 j T 1 2I 2 Z p K 2 m u T jTj P T 

JV± ~~ 77 727 T 1 — 7 : — Q1 r o tt 727 T x - [° z ) 

nibm [It + {Auj 7w ra ) J ^7 fem [7f. + Auj 2 \ [7^ + (Acj =f u m ) \ 

The =Fl term in the back reaction noise parameter depends on whether the filter is centered 
at uj s = Up + uj m or uj s = u p — u m and corresponds respectively to 'phase preserving' or 
'phase conjugating' detection as discussed in Caves.— In the limit Iq — > and or KT m — > 0, 
we see from Eqs. (179]) . (IHT]) . and (182]) that the back reaction noise amounts to doubling the 
oscillator quantum zero-point motion signal in the phase conjugating case, while the back 
reaction noise exactly cancels the quantum zero-point motion signal in the phase preserving 
case. In both cases, the noise coincides with the lower quantum bound (175]) . However, in this 
small drive/coupling limit, we do not have a detector or amplifier but rather an attenuator, 
which is of only academic interest to us. 
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Comparing the detector response (179]) and back reaction part of Eq. (1HT1) . we see that 
the mechanical oscillator behaves in the steady state as if in contact with a thermal 
bath.-'^^^'^^i'^ The back reaction of the detector on the mechanical oscillator is effec- 
tively that of a thermal bath with damping rate 7back = 1bm{R-y ~ 1) an d effective thermal 
average occupation number nb ac k defined as follows: 

7back(2< ack + 1) = -f bm N ± . (83) 

Thus, 



t* = (Ar-l)" 1 ^±-^ (SI) 



n 

The failure to approach the lower quantum bound (1751) when N± ^> 1 then translates into 
having (2n^ ack + l)7back/76m 3> 1- Thus, to get close to the bound, we necessarily require 
7back <C 7fe m ;— the back reaction occupation number ra^ ack does not have to be small. With 
the mechanical oscillator also in thermal contact with its external bath, the net damping 
rate of the oscillator is 7 nc t = 76m + 7back = R-ylbm and the net, effective thermal average 
occupation number n net of the oscillator is defined as follows: 

7net (2n± t + 1) = lbm $n{Rj*) m ) + 1] + 7 bac k (2^ + l) . (85) 

Thus, 



""net -"-7 



n(R w u m ) + ]- + ^N± 



\ (86) 



From Eq. (1781) . we see that depending on the detuning parameter Acu = uj v — ujt, the 
damping rate of the oscillator due to the detector back reaction can be either negative or 
positive. Specifically, positive damping requires the following condition on the detuning 
parameter: 

Acj < \ C \ 2u T K Tm _ 2I l Z P K Tm^TlpT ^ 

ixoj m huj m (7y + Au 2 ) 

B. Displacement sensitivity 

In the absence of a classical force acting on the mechanical oscillator, from Eq. (179"1) the 
mechanical oscillator thermal noise displacement signal spectral density takes the familiar 
Lorentzian form: 

= 7 lB \/ B ^2^T [2n {R ^ m) + 1] ' (88) 
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In order to be able to resolve this mechanical signal, the detector noise (|8T]1 referred to the 
mechanical oscillator input must be smaller than flHHl) . The detector noise spectral density 
at the input is 



J^ 7 T&m 



\c\ 2K TmUTlT 



T^lbm [lT + ( Aw =F U m ) ] 



^ 2nR 7 [ lT + (Au±u m ) 2 ] \ h 

\c\ 2 K Tm u T >yp T J 2mR UJ uj m ' 

where the first term on the right hand side is the back reaction noise acting on the mechanical 
oscillator and the second term is the output, probe line zero-point noise referred to the input. 
Note that the noise has been evaluated at ui = Up^R^Um, the maximum of the back reaction 
Lorentzian. 

If the detector output is to depend linearly on the mechanical oscillator signal input 
(i.e., function as a linear amplifier), then back reaction effects must be small. In particular, 
we require that 7b ac k <C Jbm, i-e., J? 7 ~ 1. With |c| being proportional Jo, we see from 
Eq. (I89I) that increasing the drive current amplitude J increases the back reaction noise, 
but decreases the probe line noise referred to the mechanical oscillator input. Thus, there is 
an optimum Jo such that the sum S x \ „ ois e is a minimum. Making the approximation J? 7 = 1 
and Ru = 1 in Eq. fl89l) and optimizing with respect to |c|, we find 



h 

S x (u = uj p ± RuU m ) | noise _ optimuin = 



OJr. 



\2 



OJr. 



\2 



■ (90) 



From Eq. ( |90l) . we see that the noise is further reduced if (i) the dominant source of trans- 
mission line mode dissipation is due to energy loss through the coupled probe (information 



gathering) line:— 7t ~ 7 p t; (ii) the detuning frequency is chosen to be Aco = TyTr^" u mi 
where the minus (plus) sign corresponds to phase preserving (conjugating) detection. With 
this detuning choice, the condition i? 7 w 1 requires (tu m /7T) 2 <C 1 and so the minimum 
detector noise is 

S x (u = u p ± R^ m ) | noisc „ optimum = — ^— [2THO ((^ m /7T) 2 )] , (91) 

mu; m 7 bm 



where in order to determine the O ((a; m /7r) 2 ) term, the full form of i? 7 given in Eq. ( ITSil 
must be used in Eq. (1891) when optimizing. Comparing with Eq. (188j) for the signal noise, 
we see that to leading order the detector noise effectively doubles the zero-point signal in 
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FIG. 2: Displacement detector noise spectral density (solid line) and lower bound (dashed line) 
versus drive current amplitude. The noise densities are evaluated at u = ujp + R^uJm, corresponding 
to phase preserving detection. 

the phase preserving case. This exceeds the lower bound on the detector noise derived from 
Eq. (I75p . which is zero to leading order in the phase preserving case. 

We now numerically evaluate Eq. (I89j) for the detector noise. The feasible example 
parameter values we use are:— B ext = 0.005 Tesla, Z p = 50 Ohms, u T /2n = 3x 10 9 s -1 , Qt = 
u; T /(27 T ) = 100, 7t = 9.4 x 10 7 s -1 , Z osc = 5 fim, A = 1 (geometrical correction factor), m = 
10~ 16 kg, uj m = 2.5 x 10 7 s _1 , and Qb m = ^m/i^lbm) = 10 3 . These values give a mechanical 
oscillator zero-point uncertainty Ax zp = 1.45 x 10~ 13 m, a zero-point displacement noise 
hf (mUmjbm) — 3.4 x 10 -30 m 2 /Hz, and a dimensionless coupling strength KT m = —1.1 x 
10 -5 , where we assume that in the expression ( 1221) for KTm, $ext can be chosen such that 
the dimensionless factor 4n ^° n tan (7r$ ext /<l>o) sec (7r$ e xt/$o) ~ 1 (matching condition). We 
also suppose that 7t ~ 7 p t, i.e., the transmission line mode damping is largely due to the 
probe line coupling. 

Fig. [2] shows S x (u = uj p + i?o;^m)| noise x mujm'ybm/fr and also the lower bound on the 
detector noise that follows from Eq. (1751) for phase preserving detection. Note that the 
minimum detector noise is approximately 0.8 h/(muj rn 'ybm)- Thus, for this example, the 
next-to-leading O {{oj m /^T) 2 ) term in Eq. (19T1) is approximately —0.2. Note also that the 
detector noise coincides with the lower bound in the small drive limit. 
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C. Force sensitivity 



Consider a monochromatic classical driving force with frequency ujq ~ R w uJ m acting on 
the oscillator: F ext (ui) = F 5(u — lj ). The force signal spectral density is then Sp(uj)\ signal = 
Fq5{uo — u)q). For force detection operation, the mechanical oscillator is included as part 
of the detector degrees of freedom. From Eqs. (I79ti82p . the force noise spectral density 
evaluated at to = uj p ± uj q is 



| \2j^2 ^j2 ,y 

S F {uj = Up± u ) | noise = 2mhw w n hm { 2n (u Q ) + 1 =F 1 + 




\2 

,00, 



^ bm \c\ 2 K^ m u^ pT 




(92) 



Comparing the displacement noise fl89|) with the force noise fl92l) . we see that the latter 
includes the additional 2m/ia; m 7b m [2n (uj ) + 1] mechanical quantum thermal displacement 
noise term. Since the mechanical oscillator forms part of the force detector, it need not 
necessarily be weakly driven and or weakly coupled to the transmission line; as explained in 
Sec. IHI A| the present analysis employs a linear response approximation for force detection, 
not displacement detection. Thus, in determining the optimum I (and or K^rn) and Auj 
such that SVInoise i s a minimum, we should not assume a priori the restrictions R y , R w « 1. 




lo(10-° A) 



FIG. 3: Force detector noise spectral density versus drive current amplitude for detuning Auj = 
(solid line), Auj = —5uj m (dashed line), and Au = — Wu) m (dotted line) . The noise densities are 
evaluated at ui = ui p + R^uJm, corresponding to phase preserving detection. 
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Fig. [3] shows the results of numerically evaluating the force noise spectral density given 
by Eq. ([9*2"]) for phase preserving detection (a? = u p + u ) and a range of detuning values. 
The same example parameters are used as in the above displacement sensitivity analysis, 
with n(u Q ) = and u = R u oJ m . The force noise is expressed in units 2m}TUJ rn ji )m = 6.6 x 
10~ 39 N 2 /Hz. Note that the minimum force noise is exactly 2 in these units, independently 
of the detuning, with the minimum occuring at larger Jo values as the detuning is made 
progressively more negative. 



D. Back reaction cooling 

From Eq. (|86p . we see that the net, thermal average occupation number n net of the 
mechanical oscillator's fundamental mode decreases as i? 7 increases. Thus, by increasing 
the drive and or coupling strength such that 7b ac k 3> Ibm, the mechanical oscillator can be 
effectively cooled at the expense of increasing its damping rate.-" 



Consider sufficiently negative detuning such that — Au 3> \c\ 2 u^K^ m / (iiu m ) [see Eq. (JHj 
Substituting definition f|751) for i? 7 and definition f[8~2"j) for N + into Eq. (1561) and supposing i? 7 
is large enough that we can neglect the external damping term 7^ m , we obtain approximately 
for the phase preserving case: 

n nct ~ D + n back; l yd J 

itry 

where 

_ 1 2 T +(Au + u m ) 2 _ 1 
nback ^ 4Ao^ 2' 194 j 



This expression agrees with that derived in Ref . [26|, apart from the 1/2 which is simply due 
to a small difference in the way we define n^ ack in Eq. (!83j) . Choosing optimum detuning 



Alu = — \/7y + cu^ to minimize n^" ack in Eq. (|9~~J~) . we therefore have 

How much cooling can be achieved depends on (i) how large _R 7 can be, subject to the above 
inequality on — Am; (ii) making the ratio 7^/u; m as small as possible. 26 

Using the same example parameter values as above, but taking instead a larger but still 
realistic quality factor Qh m = 10 4 for the mechanical oscillator," the resulting numerically 
evaluated effective occupation number n+ t [Eq. fl86|) ] is given in Fig.H]for a range of external 
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FIG. 4: Net effective average occupation number of the mechanical oscillator versus drive current. 
The solid curve is for external bath temperature T = 100 mK [n(R w uj) = 523], the dashed curve is 
for T = 10 mK [n(K w w) = 52], and the dotted curve is for T = 1 mK [n(R w tu) = 4.8]. 

bath occupation numbers n(u; m ). Thus, even for small coupling strengths Kx m and drive 
current amplitudes Jo, significant cooling of the mechanical oscillator can be achieved. This 
is in part a consequence of the fact that the quality factor Qh m of the mechanical oscillator 
when decoupled from the detector is very large. 

V. CONCLUDING REMARKS 

In the present paper, we have attempted to give a reasonably comprehensive analysis 
of the quantum-limited detection sensitivity of a DC SQUID for drive currents well below 
the Josephson junction critical current I c . In this regime, the SQUID functions effectively 
as a mechanical position-dependent inductance element to a good approximation and the 
resulting closed system Hamiltonian (I24|) takes the same form as that for several other 
types of coupled mechanical resonator-detector resonator systems. Thus, the key derived 
expressions ( l69l) and ( J70l) for the detector response and detector noise are of more general 
application. 

The main approximation made in analyzing the position and force detection sensitivity, 
as well as back reaction cooling, was to limit the drive current and or coupling strength 
according to Eq. fl76l) . This allowed us to find much simpler, analytical approximations to 
the key expressions, in particular Eqs. (1791) and (I8T1) . The regime of larger drive currents and 
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or coupling strengths which exceed the limit (1761) remains to be explored. However, with 
the SQUID in mind, it is more appropriate to consider larger drive currents in the context 
of including the non-linear I/I c corrections to the SQUID effective inductance. This will be 
the subject of a forthcoming paper.— 
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